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A thermodynamic approach to derive the liquid-glass transition line in the reduced temperature 
vs reduced density plane for a monatomic Lennard-Jones fluid is presented. The approach makes 
use of a recent reformulation of the classical perturbation theory of liquids [M. Robles and M. Lopez 
de Haro, Phys. Chem. Chem. Phys. 3, 5528 (2001)] which is at grips with a rational function 
approximation for the Laplace transform of the radial distribution function of the hard-sphere fluid. 
The only input required is an equation of state for the hard-sphere system. Within the Mansoori- 
Canfield/Rasaiah-Stell variational perturbation theory, two choices for such an equation of state, 
leading to a glass transition for the hard-sphere fluid, are considered. Good agreement with the 
liquid-glass transition line derived from recent molecular dynamic simulations [Di Leonardo et al., 
Phys. Rev. Lett. 84, 6054(2000)] is obtained. 
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The liquid-glass transition in Lennard-Jones (LJ) flu- 
ids has been a subject of interest for over twenty five 
years Recently, Di Leonardo et al. Q using molecu- 
lar dynamic simulations determined the liquid-glass tran- 
sition line of a monatomic LJ system in the reduced 
density (p*) vs. reduced temperature (T*) plane. Here 
p* = pa 3 and T* = ksT/e, with p the number den- 
sity ks the Boltzmann constant, T the temperature, and 
a and e being the usual parameters of the LJ potential 
(<l>Lj(r) = 4e(cr 12 /r 12 — er 6 /r 6 ), where r is the distance). 
Based on these results, they propose an off-equilibrium 
criterion to define the glass transition temperature T g . 
This effective temperature compares rather well with the 
T g obtained from equilibrium calculations. 

In a completely different context, a long time ago Hud- 
son and Andersen S addressed the nature of the glass 
transition in monatomic liquids through an equilibrium 
calculation. In particular, in the case of the LJ fluid, they 
used the Weeks-Chandler- Andersen (WCA) perturbation 
theory of liquids and two independent indications that 
a change in properties similar to a glass transition hap- 
pens in the hard-sphere (HS) fluid at a packing fraction 
of 0.533 ±0.014. Given the then scarce amount of simula- 
tion data available to compare with, their conclusion was 
that the use of the HS fluid as a reference fluid within the 
WCA scheme was adequate to derive the glass-transition 
line of monatomic fluids. As they already pointed out, 
for their approach to be useful in conection with glass 
formation in real systems, it is crucial that the HS fluid 
itself undergoes a glass transition. As it is well known, 
the existence of a glass transition in the HS fluid has been 
a debatable issue for a long time, but evidence from var- 
ious sources seems to suggest that this is indeed the case 

ill 1 1 p. 

The main aim of this letter is to assess whether an equi- 
librium calculation, such as the one carried out by Hud- 
son and Andersen H, can still be useful to describe the 



new data on the liquid-glass transition line of Lennard- 
Jones fluids. In a recent paper p!l| , we have reformu- 
lated the most popular schemes of the perturbation the- 
ory of liquids, namely the Barker-Henderson |12j, the 
variational Mansoori-Canfiled/Rasaiah-Stcll |i~3|, |14| and 
the WCA schemes, using the HS fluid as the refer- 
ence fluid. Our study focussed on the equilibrium prop- 
erties (equation of state, critical point, liquid and solid 
branches of the reduced temperature vs. reduced density 
curve at coexistence and radial distribution function) of 
the Lennard-Jones fluid. All the calculations related to 
the different perturbative schemes that we reported in 
ref. fli"|| rely on a known analytical rational function 
approximation (RFA) to the radial distribution function 
(rdf) gHs( r ) for the HS fluid developed by Yuste and San- 
tos p5[ . Their main idea is to propose that the Laplace 
transform G(t) = C\rgns(r)\, where t is the Laplace 
transform variable, may be approximated by 
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where r\ — w/6pd 3 is the packing fraction, d being the 
hard-sphere diameter and $(i) is a rational function of 
the form 

_ 1 + Sit + s 2 t 2 + S 3 t 3 + S 4 t A 

The coefficients Si and Li are algebraic functions of r\. 
They are determined by imposing two physical restric- 
tions to the rdf ]iq| , namely 

1. The first integral moment of hns(f) = 9hs{t) — 1 
is well defined and non-zero. 

2. The second integral moment of hnsir) must guar- 
antee the thermodynamic consistency of the com- 
pressibility factor Zhs — vI^P^bT) ( p being the 
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pressure) and the isothermal susceptibility \hs 
(d(r ) Z HS )/dr,)- 1 . 

These two conditions readily imply |l5|, that 
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H-e, Z PY = ^1 and Xpy = gz^ 
compressibility factor and isothermal susceptiblity aris- 
ing in the Percus-Yevick theory. In order to close the 
approximation, a given equation of state for a HS fluid 
i.e. an explicit expression for the compressibility factor 
Zhs is needed. 

All the perturbation schemes mentioned above intro- 
duce an effective (in general density and temperature de- 
pendent) diameter of the spheres as a fitting parameter 
to adjust some of the thermodynamic properties of the 
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denote the 



system of interest with respect to those of the reference 
HS system. Once this effective diameter is determined, 
one can infer by inversion the values of the temperature 
and density in the real system that correspond to a given 
packing fraction of the HS fluid. This fact was used in 
rcf. |ll| to determine the liquid and solid branches of the 
reduced temperature vs. reduced density curve at coex- 
istence for a LJ fluid from the knowledge of the packing 
fractions for the fluid-solid transition r]F-s = 0.494 and 
the solid-fluid transition t]s-f = 0.54 in the HS fluid, 
respectively |T?j ]. In a similar fashion, the liquid-glass 
transition line for the LJ system in the p* — T* plane 
may be derived from the simple relationship 



-p*d* a (p*,T*)= Vg , 
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where ?7 S is the packing fraction corresponding to the 
glass transition in the HS fluid and we have introduced 
the reduced diameter (in units of <r)denoted by d*. 

In order to proceed with such derivation, one has to 
specify the perturbation scheme and to know the value 
of 77 9 . As already stated above, Hudson and Andersen 
used the WCA and took 77^ = 0.533 ± 0.014. In our case, 
we will consider two different equations of state for the 
HS fluid that allow us to calculate in a self-consistent way 
the value of 77 g , and take the MC/RS scheme. 

The first such equation of state is the Pade(4,3) con- 
structed from the knowledge of the first eight virial coeffi- 
cients jl8|, |l^] . The compressibility factor corresponding 
to this equation of state is given by 
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Note that, as was pointed out in previous work[g[ [20| . 
within the RFA the Pade(43) leads to a glass transition 
in the HS fluid at ?7 g = 0.5604. It should also be borne 
in mind that this equation of state has a simple pole 
at a packing fraction very close to the fee close-packing 
fraction [p|. 

On the other hand, the second equation of state is an 
ad-hoc approximation constructed in the following way. 
Since the well known Carnahan-Starling(CS) pl[ equa- 
tion of state has been shown to be very accurate through- 
out the complete fluid region and even in a small density 
range within the metastable regime but has a pole at 
7/ = 1 which is clearly unphysical, we demand that the 
new compressiblity factor Z prop has a pole at the ran- 
dom close packing fraction r] rcp (up to this point taken 
as a parameter) and that the first eight coefficients in its 



series expaexpansionnsion coincide with the same coeffi- 
cients in the series expansion of the Carnahan-Starling 



compressibility factor, 
the form 



Thus, Zp rop (r)) is taken to be of 
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where the coefficients dj (i — 1,...,5) and bj (j = 1,2) 
depend on the value of ?7 rC p. In order to determine this 
value, we further assume that, according to the crite- 
rion that has been used earlier^, |2(|, a glass transition 
in the HS fluid (considered to be a second order phase 
transition) occurs at 77 g (also unspecified at this stage) 
if 54 vanishes for r] g . The above conditions lead to the 
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following set of equations 
= 0, 
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where A is a constant and in writing the compressibility 
factor for the glass we have taken the form proposed by 
Speedy §, §. The set of equations @-(||) together 
with (|ll[) allow us to determine rj rcp , A and r\ g in a self- 
consistent way. The results are rj rcp = 0.6504, A = 2.780 
and ri g = 0.5684, which are well within the range of pub- 
lished values for these quantities. In turn the value of 
rjrcp leads to the explicit form of Z prop (r)), namely 
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In the range < rj < 1 this equation of state only presents 
a simple pole at r\ = r\ rcp . 

In order to illustrate the numerical accuracy of these 
two equations of state in the metastable region, in Fig. 
[I] we display the inverse of the contact values of the rdf 
(g(d + )^ 1 = At]/ (Zhs — 1)) derived from them as a func- 
tion of the packing fraction and compare the results with 
the simulation data obtained by Rintoul and Torquato 
|p2| . In this figure we have also included the predic- 
tions of the corresponding equations of state for the glass 
adopting the form suggested by Speedy H Q]. 

For the sake of making the paper self contained, we 
now recall the condition required to derive the tempera- 
ture and density dependent diameter within the MC/RS 
perturbation scheme (for details see refs. [O, [l3| Q). 
Here the effective diameter is chosen to minimize the 
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FIG. 1: Inverse of the contact values of the rdf as derived from 
Z43 (continuous line) and Z prop (dashed line) as a function of 
77. The region corresponding to the glass has been obtained 
from Z „ias s and in the case of Z43, A = 2.765 and Tj rcp = 
0.6448 fed]. The full symbols are simulation results of ref. 
and the arrows indicate the location of r\ g for each case. 



Helmholtz free energy of the perturbed system. To first 
order in /3 = l/T* this diameter is determined from the 
equation 111) 
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Clearly in solving this equation it should be understood 
that Zhs and <&(£) have been expressed as functions of 
p* and d* . 

Using the effective diameters determined from eq.(|l6|) 
with either Z43 or Z prop and the condition given by eq. 
(0) (with the rig value corresponding to each equation 
of state) we have determined the liquid-glass transition 
lines in the p*-T* plane for the LJ fluid. These are shown 
in Fig. where we have also included the recent simu- 
lation results of Di Leonardo et al. . As clearly seen 
in the figure, not only the qualitative trend observed in 
the simulations is reproduced, but also the quantitative 
agreement is rather good, specially for the case of Z prop . 
A remarkable aspect of these results is that they have 
been derived self-consistcntly with no free parameters. 

One question that immediately arises is to what extent 
the above results depend on the MC/RS perturbation 
scheme and the choices for Zhs- Concerning the first 
issue, which is particularly relevant in view of the fact 
that the previous calculation by Hudson and Andersen 
H used the WCA, we have checked that the performance 
of this latter scheme is much poorer. This may reflect 
that the choice of minimizing the Helmholtz free energy 
of the perturbed system to derive the effective diame- 
ter may be more astringent than equating the isothermal 
compressiblities of the actual and the HS reference sys- 
tem. Something similar applies to the Barker-Henderson 
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Furthermore, the application of the same approach for 
other monatomic fluids may be laborious, but in princi- 
ple should follow the same steps. 



FIG. 2: Liquid-glass transition lines as obtained with the 
MC/RS perturbation scheme using Z43 (continuous line) and 
Zp r0 p (dashed line). The points are the simulation data of Di 
Leonardo et al. in ref Q 



perturbation theory in which the lack of density depen- 
dence of the effective diameter does not reproduce the 
desired trend. As far as the second issue is concerned, 
we have also performed calculations taking the CS equa- 
tion of state (which does not lead to a glass transition 
in the RFA formulation) and adjusted the value of the 
packing fraction rj g until good agreement with the sim- 
ulation data of Di Leonardo et al.B was obtained. In 
fact for r\ g = 0.55 we get results pretty close to those 
derived using Z prop . The above suggests that if the ex- 
act equation of state of the HS system (including a glass 
transition) were available then the use of perturbation 
theory and the RFA approach would yield a very accu- 
rate description of the liquid-glass transition line of the 
LJ fluid. 

In conclusion, in this paper we have shown that an 
equilibrium approach to the glass transition in LJ flu- 
ids, much in the same spirit as discussed by Baeyens and 
Verschelde |j) in the case of HS fluids, is wholly compat- 
ible with the molecular dynamics simulation results for 
the liquid-glass transition line reported recently |2| . Per- 
haps the key aspect of our approach is its self-consistency 
and the provision of a theory with no free parameters. 



t 

[1 

[2 
[3. 

[4; 

[5: 

ir 

[8 
[9 
[10 

[11 

[12 

[13 

[14 
[15 
[16 

[17 
[18 
[19 
[20 

[21 

[22 



mrp@cie.unam.mx 




malop ez @ servidor . unam . mx 



Electronic address: 
Electronic address: 
See for instance, A. Rahman, M. Mandell, and J. Mc- 
Tague, J. Chem. Phys. 64, 1564 (1976). 
R. Di Leonardo, L. Angelani, G. Parisi, and G. Ruocco, 
Phys. Rev. Lett. 84, 6054 (2000). 

S. Hudson and H. C. Anderssen, J. Chem. Phys. 69, 2323 
(1978). 

J. D. Weeks, D. Chandler, and H. C. Andersen, J. Chem. 
Phys. 53, 149 (1971). 

R. J. Speedy, J. Chem. Phys. 100, 6684 (1994). 

B. Baeyens and H. Verschelde, Z. Phys. B 102, 255 

(1997). 

R. J. Speedy, Mol. Phys. 95, 169 (1998). 

M. Robles, M. Lopez de Haro, A. Santos, and 

S. Bravo Yuste, J. Chem. Phys. 108, 1290 (1998). 

W. van Megen and S. M. Underwood, Nature 362, 616 

(1993). 

W. van Megen, T. C. Mortensen, S. R. Williams, and 
J. Miiller, Phys. Rev. E. 58, 6073 (1998). 
M. Robles and M. Lopez de Haro, Phys. Chem. Chem. 
Phys. 3, 5528 (2001). 

J. A. Barker and D. Henderson, J. Chem. Phys. 47, 2856 
(1967). 

G. A. Mansoori and F. B. Canfield, J. Chem. Phys. 51, 
4958 (1969). 

J. Rasaiah and G. Stell, Mol. Phys. 18, 249 (1970). 
S. B. Yuste and A. Santos, Phys. Rev. A 43, 5418 (1991). 
M. Robles and M. Lopez de Haro, J. Chem. Phys. 12, 
4648 (1997). 

J. P. Hansen and L. Verlet, Phys. Rev. 184, 151 (1969). 
E. J. J. van Rensburg, J. Phys! A 26, 4805 (1993). 
I. C. Sanchez, J. Chem. Phys. 10, 7003 (1994). 
S. Bravo Yuste, M. Lopez de Haro, and A. Santos, Phys. 
Rev. E 53, 4820 (1996). 

N. F. Carnahan and K. E. Starling, J. Chem. Phys. 51, 
635 (1969). 

M. D. Rintoul and S. Torquato, J. Chem. Phys. 105, 
9258 (1997). 



